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We solve the Navier-Stokes equations with two simultaneous forcings. One forcing is applied at a given 
large-scale and it injects energy. The other forcing is applied at all scales belonging to the inertial range and it 
injects helicity. In this way we can vary the degree of turbulence helicity from non helical to maximally helical. 
We find that increasing the rate of helicity injection does not change the energy flux. On the other hand the 
level of total energy is strongly increased and the energy spectrum gets steeper. The energy spectrum spans 
from a Kolmogorov scaling law for a non-helical turbulence, to a non-Kolmogorov scaling law for 
a maximally helical turbulence. In the later case we find that the characteristic time of the turbulence is not the 
turnover time but a time based on the helicity injection rate. 

We also analyse the results in terms of helical modes decomposition. For a maximally helical turbulence one 
type of helical mode is found to be much more energetic than the other one, by several orders of magnitude. 
The energy cascade of the most energetic type of helical mode results from the sum of two fluxes. One flux 
is negative and can be understood in terms of a decimated model. This negative flux is however not sufficient 
to lead an inverse energy cascade. Indeed the other flux involving the least energetic type of helical mode is 
positive and the largest. The least energetic type of helical mode is then essential and cannot be neglected. 

PACS numbers: 47.27.-i, 47.27.Ak, 47.27.Gs 


Considering the case of three-dimensional homogeneous 
and isotropic turbulence, Kolmogorov HI] assumed the exis¬ 
tence of a range of scales, the so-called inertial range, in which 
the viscous dissipation can be neglected. In absence of dissi¬ 
pation the kinetic energy is a conserved quantity. In spectral 
space the flux of kinetic energy is constant, leading to an en¬ 
ergy cascade from large to small scales, provided energy is 
injected at large scale. From a straightforward dimensional 
analysis the spectral density of kinetic energy E{k) can be ex¬ 
pressed in terms of the energy injection rate e and wave num¬ 
ber k, 

E(fc) a £2/3fc-5/3. (1) 

In absence of viscosity not only energy is a conserved quan¬ 
tity but also helicity ytot), which is defined as 

H(f) = f u(x, f) • V X u(x, t) dV, (2) 

Jv 

where u(x, t) is the velocity field at position x and time f, and 
integration is made over the volume V. Similarly to energy, 
the helicity conservation is equivalent having a constant he¬ 
licity flux. Provided energy and helicity are both injected at 
large scale, helicity is expected to cascade jointly with energy 
in the inertial range, obeying the following scaling law IJ-l^ 

H{k) oc ;;e-V3fc-5/3 ( 3 ) 

where H{k) is the spectral density of helicity and t] the injec¬ 
tion rate of helicity. Contrary to enstrophy in two-dimensional 
turbulence, helicity is not sign-defined and therefore not re¬ 
puted for having any influence on the energy spectrum, letting 
the scaling law ([TJ unchanged ||3]. 


Though the simultaneous scaling laws ([B and Q are char¬ 
acteristic of the so-called helical turbulence |@], they cannot 
be justified from dimensional grounds like Kolmogorov did 
for the non helical turbulence. Indeed the problem now con¬ 
sists in five variables E{k),H{k),k,e and t] and only two di¬ 
mensions, length and time. Applying the Fl-theorem |0] and 
assuming that E{k) and H{k) obey some scaling laws, we And 

E(k) oc £7/3-«^«-5/3^-fl^ ^ ^4/3-fc^fc-2/3^-b ( 4 ^ 

where a and b are two free parameters. Therefore we need 
additional constraints to derive the power laws for E{k) and 
H{k). 

One way to argue for the simultaneous k scaling laws 
(HJ and (IB is to assume that the fluxes of energy and helicity 
nE(fc) and Tli{{k) are constant in the inertial range, such that 
Tl£{k) = £ and Tlj^{k) = rj. In addition we have to set that 
the characteristic times Tg and Th for the energy and helicity 
transfers are given by the turbulence turnover time T£ = Th oc 
(ek^)~^/^. Then estimating both energy and helicity fluxes as 

ill [11 

UEik) = kE{k)/TE{k), Tlnik) = kHik)/TH{k), (5) 

leads to ([B and (|B. In the notations of (IBi this would corre¬ 
spond to fl = b = 5/3. 

Instead we could think of spectral laws independent of £, 
leading to Jl] 

E{k) oc /;2/3fc-7/3^ H{k) a ;;2/3fc-4/3_ ( 5 ) 

In the notations of (HJi this would correspond to a = b - 1 -1 = 
7/3. Providing evidence of such scaling laws (IB, is still a 
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challenging issue and has never been observed so far in di¬ 
rect numerical simulations. Recently, a step forward has been 
made by solving the so-called decimated Navier-Stokes (NS) 
equations lll2[ll3n . It consists in splitting each Fourier modes 
of the velocity field in positive and negative helical modes, 
and in solving the NS equations keeping only one type of 
mode. By construction the resulting turbulence is then ex¬ 
actly maximally helical, i.e. \H{k)\= kE(k) . In such a deci¬ 
mated model helicity is still a conserved quantity, but now it 
gets the property to be sign-definite. It then plays a role sim¬ 
ilar to enstrophy in 2D turbulence, leading to a k~^l^ inverse 
cascade of energy at scales larger than the forcing scale 11211. 
In the inertial range the helicity cascade is direct, with an en¬ 
ergy scaling law E{k) oc 0. In a recent experimental 

study 0 two scaling laws have also been found, the authors 
arguing for the existence of two such opposite cascades, but 
with dominant non-local transfers leading to E{k) a k^^ and 
E{k) a fc-2. 


Here we present another strategy that does not assume any 
simplification of the NS equations and which is based on the 
fact that the scaling laws (|6]l do not depend on e. Such a e 
independence is expected as soon as < T£, then speaking 
of a helicity driven turbulence. From Q and applying the ex¬ 
act constraint |H(A:)|< kE{k), a sufficient condition for having 
Th < Te is given by 


Unik) > kllEik). (7) 

One way to satisfy such a flux condition (|7]) is to inject energy 
at large scale such that He (A:) = e and helicity at all scales 
such that nEf(A:) > ke. This is the strategy that is followed in 
the present letter using direct numerical simulation of the NS 
equations. A similar strategy has been used in 00 using 
a helical shell model. 

Using a pseudo-spectral code we solve the NS equations 

dfU = —(u ■ V)u — Vp-t vV^u-t f (8) 

where v, p and f are respectively the viscosity, the pressure 
and the flow forcing. The forcing is divided in two parts f = 
f£ ^ where is the energy forcing applied at some given 
large scale Ap^, and is the helicity forcing applied at all 
scales within the inertial range. 

Both parts of the forcing and are delta-correlated in 
time and divergence free. Following 0 they are defined 
such that the power input comes from the force-force cor¬ 
relation only and not from the velocity-force correlation. In 
spectral space this corresponds to 

<-ff + c.c. = 0 (9) 

Up -ff-tc.c. = 0 (10) 

where fp and are the Fourier coefficients of and f^. 

For we use the exact same forcing as in 0 with a force- 
force correlation given by 

|ff|2 = F{k)/2nk\ (11) 


where F(k) obeys to a Gaussian distribution around k = fcp. 
As in MITT F(k) is defined as inversely proportional to the time 
step of the computation, in order to guarantee an injection rate 
of energy which is independent from the value of the time 
step. The level of helicity injected by is not controlled a 
priori, but the results show that it is statistically insignificant. 

In order to inject helicity the forcing has to satisfy, in 
spectral space, 

(V X u)p ■ -I- Up ■ (V X f^)p -H c.c. = ri{k), (12) 

where t]{k) is a helicity injection rate per unit volume. We 
take 


?7(k) = 0 

?/(k)= T]o{\k\/kp) 


for |k|< kp 
for |k|> kp 


(13) 

(14) 


with a = 2.2 in order to have a spectral density of helicity 
injection rate |kp;/(k) almost flat. Of course such a forcing 
extending on the whole inertial range might change the inter- 
mittency properties of the turbulence llSll . However we And 
that the level of dissipation with and without is statistically 
unchanged. Finally two issues have to be clarified, both re¬ 
lated to the fact that the energy power coming from the force- 
force correlation of is not controlled a priori and that we 
need to keep it at a level sufficiently lower than the one in¬ 
jected by f^. These rather technical issues are detailed in the 
Appendix. 

Applying a classic criterion 0 in order to ensure the reso¬ 
lution of a sufficiently large range of dissipation scales, taking 
a grid of 256^ points and setting v = 2 10“^ and = 100, 
where R\ is the Reynolds number based on the Taylor mi¬ 
croscale, leads to an energy injection rate e w 0.2 and a forc¬ 
ing wave number kp w 2.2. Finally all subsequent results 
correspond to statistically steady states. 

In Figure [U top and bottom, the spectral density of energy 
E{k) and relative helicity H{k)/{kE{k)) are represented for 
five values of helicity injection t]Q, ranging from non-helical 
turbulence (t/q = 0) to maximally helical turbulence {tjq = 5). 
Clearly, increasing f/o steepens the energy spectral density at 
large scales, with a scaling law varying from k~^l^ for non¬ 
helical, to for maximally helical turbulence (top of Fig- 
ure[I]l. For t]Q =/ 0 a well defined spectrum of relative helicity 
is obtained with a rather flat part. For tjq = 1 and t]Q = 5 
the relative helicity is about unity over an extended range of 
scales, showing that the turbulence is close to a maximally 
helical state (bottom of Figure[T]i. 

The fluxes of energy and helicity, nE(A) and nEf(A), are 
plotted in top and bottom of Figured for again the same five 
values of j/q- In the top figure we see that He (A) is almost 
independent of tjq, showing that the spurious energy injection 
produced by the helical forcing is small compared to the 
energy injected by f^. On the other hand in the bottom fig¬ 
ure we see that nEf(A) is getting higher when increasing the 
value of //Q. Therefore we conclude a posteriori that the injec¬ 
tions of energy and helicity are well prescribed by and 
respectively. 
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FIG. 1: (Color online) Spectral density of energy (top) and relative 
Felicity (bottom), for five values of the Felicity injection rate: rjQ = 
0,0.05,0.3,1 and 5. In the top figure the energy is normalized by 
and the red dashed curve corresponds to 


Relying on (|5]l and knowing the flux and spectral density 
of energy and Felicity, we can calculate the two characteristic 
times, TE(k) and Th(A:), in order to determine which one is the 
smallest and therefore which one controls the turbulence. In 
Figure[3]the ratio T]^(k) /Te(A:) is plotted for j/g = 0.05,0.3,1 
and 5. For sufficiently large values of t]Q, typically J/g = 1 and 
rjQ = 5, we see that for k G [3,12] TEf(fc)/T£(fc) < 1, sug¬ 
gesting a turbulence governed by the Felicity injection rate. 
On the other hand for low values of tjq, typically tjq = 0.05 
and = 0.3, and in the same range of scales we find 
'^h(^)/'1e(A:) > 1 , suggesting a turbulence governed by the 
energy injection. 


Up to now it has been demonstrated that injecting a suffi¬ 
ciently high rate of Felicity over the whole inertial range of a 
turbulent flow leads to a scaling law for the energy spec¬ 
tral density and that the characteristic time of such maximally 
helical turbulence is the one based on the Felicity injection 
rate. Apart from studies assuming drastic simplification of the 
NS equations \12, 13,15, 3, this is the first direct numerical 
simulation in which Felicity is shown to have some effect on 
three-dimensional homogeneous and isotropic turbulence. 

Now the question arises how our results fit in with the sce¬ 
nario described by the decimated model lll2l I Of . As the in¬ 
jection of positive Felicity is made at all scales and at each 
time step, we expect a strong dominance of the positive he- 


FIG. 2: (Color online) Flux of energy He (top) and Felicity 11^ 
(bottom) for the same five values of tjQ as in Figure[T]and same color 
code. 



FIG. 3: (Color online) Ratio T]^{k)/TE{k) versus k for rjo = 
0.05,0.3,1 and 5 and same color code as in Figure[T] 


lical modes compared to the negative helical modes. Then 
according to lEQ we could expect an inverse cascade of 
energy. However the energy fluxes plotted in Figured are al¬ 
ways positive, demonstrating a direct cascade of energy. To 
clarify this paradox we now analyse our results in terms of 
helical modes decomposition. 

In Fourier space the velocity field is split into two helical 
modes per wave vector 

u(k) = u+(k)-i-u^(k) (15) 

= (k)h’^ (k) I- M “ (k)h“ (k) (16) 
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FIG. 4: (Color online) Energy spectra E'''{k) (solid line) and E^{k) 
(dashed line) for the same five values of tjq as in Figure [T] and same 
color code. Increasing tjq the curves spread from the center towards 
the two (red) dashed and solid straight lines corresponding to 


where ii^ are complex scalars and are the eigenvectors 
of the curl operator satisfying ik x = ±|k|h^ 12^ 2 lll . 
The energy spectral density of each helical mode, defined by 
E^(k) = \u^{k)f', is plotted in FigurelH for again the same 
five values of //o as in Figure [T] For each value of we ob¬ 
serve that both spectra E^{k) obey to the same scaling laws, 
again varying from for tjq = 0 to for tjq = 5. 

This is consistent with the results of Figure [T] and the relation 
E{k) = E'''(k) + E~{k). For tjQ = 0 both spectra E^{k) are 
identical as expected for a non helical turbulence. Increasing 
rjQ both spectra separate apart from each other, E'''{k) prevail¬ 
ing over E~( k) by two orders of magnitude for fjQ = 5. 

Following 122112311 we now analyse the fluxes between the 
helical modes. We denote by Tl‘^'^{k) the energy flux from the 
inside region of a u®-sphere of radius k to all wave numbers 
of u^, where a,b = zk. It is defined as 



FIG. 5: (Color online) Top: Various energy fluxes in helical turbu¬ 
lence. n^''(fc) denotes the energy flux from the inside region of a 
u^-sphere of radius k to ah wave numbers of u^, where a,b = ±. 
Bottom: The dark (light) curves correspond tofjQ = 5 ()/o = 0). The 
solid curves correspond to n+<(fc) and nz‘^(fc), the dashed curves 
to nt.^(fc) and nT^(fc). For rjg = 5: (a) corresponds to !!;[;<(fc), (b) 
to ni^(fc), (c) to nT'^(fc) and (d) to nT^(fc). 


Finally we push one step further by splitting the flux 
n;[;>(fc) into two parts n;[;>(fc) = -i- “n;[;>(fc) with 


nl<{k)= [ u‘‘{k')-FT{{u-V)u^}{k')dk' (17) 

7|k'|<lr 


^nl<{k)= / u''(k')-FT{(u±-V)u''}(k')rfk'. (18) 

i|k'|<fc 


where FT{(u • V)u^} denotes the Fourier transform of the 
non linear term (u ■ V)u^. The four fluxes n^^(A:) are rep¬ 
resented schematically in Figure |5] (top). They are plotted in 
Figure |5] (bottom) for tjQ = 0 (light curves) and rjQ = 5 (dark 
curves). The sum of these four fluxes corresponds to the en¬ 
ergy flux plotted in Figure |2] (top). For t]o = 5 the fluxes (c) 
nt.^(k) and (d) 11)7 ^(k) are much larger than for j/q = 0. In 
addition they are of opposite sign corresponding to a net flux 
of energy from the positive to the negative helical modes, bal¬ 
ancing each other at small scales (large k). 

Let us now focus on the flux (a) n;[;^(A:). By definition the 
energy flux from the inside region of a u'''-sphere of radius k 
to itself is zero. This implies that n;[;^ (k) = n;j;> (k), which 
is the energy flux from the inside region of the u'''-sphere of 
radius k to the outside of that same u'''-sphere. Now the fact 
that in Figure|5](bottom) n;j;> (k) is always positive means that 
there is a direct cascade of energy. This is in contrast with the 
inverse cascade found with the decimated model of lll2l 11311 . 


In (fTSl l denotes the energy flux from the inside re¬ 

gion of a u'''-sphere of radius k to the outside of the u'''-sphere, 
with acting as a mediator on the non linear interactions. 
Both fluxes ^n;[;>(fc) are plotted in Figure |6] for again the 
same five values of rjQ as in Figure [T| The flux ‘‘‘riiji^ ik) 
is al way s negative in agreement with the arguments given in 
lll2[ll3l] for the decimated model. However the flux (k) 
is positive and always the largest in absolute value. This 
shows that even if the turbulence is strongly positively he¬ 
lical, the presence of negative helical modes is nevertheless 
essential to give the right sign of the energy fluxes. Though 
the decimated model is mathematically appealing because it 
reproduces an exact maximally helical flow, it is eventually 
singular as in practice the existence of both types of helical 
modes cannot be avoided. This result also supports the choice 
made in helical shell models ll23n in which two helical modes 
can interact only if they have opposite helicities. 

ISTerre and LEGI are part of Labex OSUG@2020 
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FIG. 6; (Color online) Energy fluxes (solid) and ^n;J;^(fc) 

(dashed) for the same five values of 1 ]q as in Figure[T]and same color 
code. Increasing tjq from 0 to 5 the solid curves at A: = 3 decrease 
and the dashed curves increase. 


to the energy rate coming from the force-force correlation of 
f^. In other words to maintain a spurious power injected by 
at a low level it is necessary to decrease the time step when 
increasing tjQ. 

Finally we apply a clipping condition in order to prevent 
any spurious energy injection coming from singular solutions 
of equations (fTOl) and (fTSI) . Indeed as the forcing is ran¬ 
dom in time we cannot prevent the value of D given in (|2^ 
to be zero and lead to singular solutions fl(k, f) and b{k,t). 
In addition we do not want to force energy nor helicity in 
the dissipative range (corresponding to scales k > ky where 
ky » Therefore at each time step the helical forc¬ 

ing is applied provided the following condition is satisfied 


D(k, t) > |k||m-|^ 


^A + B 



(24) 
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APPENDIX: THE HELICITY EORCING 

The helicity forcing is defined in its spectral form as 

ff = kxp(k,f) (19) 

with 

p(k, f) = fl(k, f)ei(k, f) + ib{k, f)e 2 (k, f), ( 20 ) 

ei(k, I) and e 2 (k, I) being two unit vectors with directions 
changing randomly at each time step. The resolution of equa¬ 
tions ([Tol l and (fT^ leads to 

=- 4D|kp (^(»A)/k,e2) (21) 

HKt) = ^^(5R(u,),k,ei), ( 22 ) 

with 

D(k,f) = (5i(ufc),k,ei) (K(ujt)- 62 ) 

+ k, 62 ) (5(ujt) • ei). (23) 

As stated above, the e nerg y forcing is inversely propor¬ 
tional to the time step lO . Conversely does not depend 
on the time step, implying that the level of energy rate which 
comes from the force-force correlation of is proportional to 
the time step. Therefore provided the time step is sufficiently 
small, the energy rate coming from the force-force correlation 
of can be maintained at a sufficiently low level compared 


where A, B and /3 are positive constants whose values depend 
on }]o. 
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